Coherence Peak in the Spin Susceptibihty from Nesting in Spin- Triplet 
Superconductors; A Probe for Line Nodes in Sr2Ru04 
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We study the dynamical spin susceptibility x{<ly^) for spin-triplet superconductivity. We show 
that a large peak at cj = 2A appears in ImXzz{Q,(^), where z is the direction of the d vector for 
triplet pairing, if Fermi surface has a nested part with the nesting vector Q and the order parameter 
is -|-A and —A in this part of the Fermi surface. If there are line nodes in the nested part of the 
Fermi surface, a peak appears in either ImXzz(Q,(^) or Imxn — (Q,aj), or both, depending on the 
perpendicular component of the nesting vector. The comparison with inelastic neutron scattering 
experiments can determine the position of the line nodes in triplet superconductor Sr2Ru04. 

PACS numbers; 74.25 Ha, 74.70 Pq, 61.12 Bt 
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I. INTRODUCTION 

The superconductivity in Sr2Ru04" has been revealed 
to be unconventional by many experiments; it is spin- 
tripleliSii, it breaks time-reversal symmetr}*4, and its en- 
ergy gap has line nodes^iSi^. 

The triplet superconductivity with horizontal line 
nodes of the energy gap has been proposed^ to ex- 
plain these experiments. The absence of the angle de- 
pendence of the thermal conductivity within the a-b 
pian o — ' J - ' ^ ' shows that the line nodes run horizontally on 
the Fermi surface. The angle resolved ultrasound attenu- 
ation experimenlii is compatible with the horizontal line 
nodes. 

The Fermi surface of Sr2Ru04 consists of three cylin- 
drical surfaces named a, /3 and ^iiSii^, as predicted by 
the band calculationi^ii^. The hybridization of the dxz 
and dyz orbits of ruthenium with the pj^ orbit of oxygen 
makes two one-dimensional bands, if the mixing of these 
bands is neglected. If small mixing is taken into account, 
we get hole-like a and electron-like /3, but the nesting of 
the Fermi surfaces survives as predicted^^ and confirmed 
by the inelastic neutron scattering experimen'feii. The 7 
surface is constructed by the dxy orbit of Ruthenium and 
is two-dimensional. 

Nomura and Yamadai^ have studied the two- 
dimensional three-band Hubbard model in the third- 
order perturbation theory and obtained that triplet 
superconductivity is stabilized mainly in the two- 
dimensional band. They also obtained the line-node- 
like power-law behavior in the temperature-dependence 
of the specific heat due to the vertical node-like energy 
gap in the a and P bands. In their treatment the momen- 
tum dependence of the gap is determined at the transi- 
tion temperature Tc. The spin-triplet state such as 



or 



d(k) = zA(cos^ kx — cos^ ky){sinkx + isinky 



(2) 



where d(k) is the d vector, should be mixed with 

z A' {sin ky + i sin kx) , (3) 



or 



zA'(sin kx — i sin ky 



(4) 



d(k) = zA sin kx sin fcj,(sin k^ + i sin ky), 



(1) 



at T < Tc and the vertical line nodes will disappear^- 

Zhitomirsky and Rice^* proposed the mechanism for 
the horizontal line nodes that while the active band has a 
full energy gap, the interband proximity effect makes the 
horizontal line nodes in the passive band. They assumed 
that the two-dimensional band is active and that the line 
nodes are in one-dimensional bands. 

On the other hand. It has been shown that spin-triplet 
superconductivity is induced in the quasi-one dimen- 
sional bands by the antiferromagnetic spin fluctuation 
if the antiferromagnetic spin fluctuation is anisotropic in 
the spin spaceSSiSi. If this is the case the one-dimensional 
bands are expected to have a full energy gap and the two- 
dimensional band has line nodes. 

Therefore, the determination of the gap structure is 
important to understand the mechanism of unconven- 
tional superconductivity in Sr2Ru04. Bulk measure- 
ments such as specific heat, NMR and thermal conduc- 
tivity cannot distinguish which part of the Fermi surface 
has the line nodes, a, (3 or 7. 

In the inelastic neutron scattering experiments, imagi- 
nary part of the dynamical spin susceptibility, Im x(q, w), 
is observed, from which we can get the information of the 
superconducting order parameter. The so-called 41 meV 
peak in the spin-singlet d-wave superconducting state of 
YBa2Cu30T^2i2i212^ has been observed. Many theo- 
retical studieo^^.^gi^9i30,3i,32,33,34,35,36,37 ^^^^ ^cen done 

to explain the peak structure in inelastic neutron scat- 
tering. 
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Dynamical susceptibility for the spin-triplet supercon- 
ductivity has been studied theoretically^Si^S. Recently, 
the resonance peak in Imx(q, oj) is shown to be a sign of 
the triplet superconductivitji^MiiiS. The order parame- 
ter assumed in these papers, however, does not seem to 
be consistent with experiments in Sr2Ru04. 

In this paper we show the general form of the dy- 
namical spin susceptibility XijCQ? ^) i'^ the unitary states 
of the triplet superconductivity and calculate it in the 
system which has the nested Fermi surface with and 
without line nodes in the nested part of the Fermi sur- 
face. We show that the position of the line nodes can 
be determined by the inelastic neutron scattering exper- 
iment, which observes the imaginary part of the dynam- 
ical susceptibilit y^'^'^'^'"^'^ 



II. DYNAMICAL SPIN SUSCEPTIBILITY OF 
SPIN-TRIPLET SUPERCONDUCTIVITY 

The dynamical spin susceptibility is given bjiS^i^MSiii 

x°i(q,«^m) = - ^ 

xT^Tr (^Q;jG(k, ie„)djG(k + q, ie„ + iujmij ,(5) 

n,k 

where LOm — 2rmrT and e,i = (2n + 1)t:T are Matsubara 
frequencies (m and n are integers), and d and G(k, ie„) 
are the 4x4 Nambu representation of the spin and Green 
function, respectively, i.e.. 



cr* 
(jya'aV 



1 + Pz , l--P±„v„z„y 
(7 -i z (J (T 



2 2 
where cr* (i = x, y, or z) is a Pauli matrix, and 



G(k,ie„) 



/G(k,ze„) F(k,ze„) 
\F\\,,ie^) -G(~k,-ze„) 



(6) 



(7) 



The 2x2 matrix Green function G(k, ie„) and the anoma- 
lous Green function i^(k, ie„) are given as the Fourier 
coefficients for 



and 



G„,/3(k, r) = -(T,ak„(T)at (0)) (8) 



F„,/3(k,T) = (T,aka(T)a_k/3(0)), (9) 



respectively. 

In this paper we consider the weak coupling theory for 
the spin-triplet superconductivity. We take account of 
the interaction U by the random phase approximation 
(RPA), 



(10) 



but the essential properties such as a peak in Imx^^ (q, uj) 
is already seen in the absence of the interaction effects. 

The order parameter is given by the d vector as 
/S.af){k) = i{{d{k)-(T)(jy)a0- We study the unitary states, 
d*(k) X d(k) = 0, in this paper, since experimental re- 
sults can be explained by the unitary states. Then we 
perform the summation over n in Eq.©, and get 



k 

+ cL'pp' (k, q)^^+^ (k, q, c.) (1 - /(i?k') - /(i?k)) I , 

(11) 

where 

'-'aa'Pp'y^^^) - 2 

'5a/3'<^o-/3aa' - Rc(A;^(k)A^,^,(kO) 

± „ , (,J-2j 



2i?k£'k' 



i?W(k,q, 



1 



1 



ii^k ± i'k' — UJ — iV' 



(13) 



k' = k + q, and the analytic continuation iujm uj + iV 
with F +0 has been done. 

We take d vector parallel to the z axis as indicated by 
experiments^'^, to obtain Xij = for i ^ j and 

x° (q,w) 

= ^ E l^i^^ (k. q)^*"^ (k, q) (/(^^k' ) - fiE^)) 

k 

+ ci-)(k,q)7^W(k,q) (1 - fiE^,) - /(E^)) |,(14) 



where 



GW(k,q) = G(±)(k,q)-Gi=^J(k,q) 

1 ^CkCk'~Re(d:(k)4(k')) 



= - ± 
2 



(15) 



and 



G(^)(k,q) . (16) 

Depending on the sign of Re(c?*(k)(iz(k')) in the co- 
herence factor, cl~^ in Eqs. ((T3|l and a peak will 

appear in either Im^n (q, i^) or hnXzzi^jUj), i.e. if 

Re((i*(k)c?^(k')) > (< 0), a coherence peak appears 
in xh (q, uj) {xzzil, w)), as we will show in the next sec- 
tion. 
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FIG. 1: Fermi surface at fc^ = and the nesting vector Q. 



III. COHERENCE PEAK FROM THE NESTED 
FERMI SURFACE 



A. case A 

First we study the case A, where we take the order 
parameters as 

di^(k) = Aisinfc^ (21) 
d22(k) = zAisinfcj,, (22) 



and 



d3,_ (k) = A3 ( sm Y cos Y + « cos y sm j 

xcos%. (23) 



In the case A, Cooper pairs are formed between the elec- 
trons on the nearest sites in the conducting plane (r^ and 
r,j-f (±a, 0, 0), ri + {0, ±a, 0)) for the one-dimensional elec- 
trons, and Cooper pairs are formed between the electrons 
on Ti and rj-|-(±a/2, ±6/2, ±c/2) for the two-dimensional 
electrons. These order parameters can be realized if the 
pairing Hamiltonian is taken as^ 



In order to study the coherence peak in Sr2Ru04, we 
take the simple three-band model, where there are two 
one-dimensional bands and a two-dimensional band, i.e., 



' £03, 



(17) 
(18) 

(19) 



-2*1 cos{k^) + eoi 

-2ti COs(fcj^) + £01 
— 2f3(cOs(fca;) -I- COs(fcj,)) 

—4^3 cos(fc2,) cos(fcy) 

where the lattice constant is taken to be 1. We set pa- 
rameters as h = 0.31 eV, eoi = -0.24 eV, ^3 = 0.44 eV, 
^3 = 0.14 eV and eo3 — —0.14 eV. The Fermi surface is 
shown in Fig. ^ We set the wave vector q to be the nest- 
ing vector Q. The z component of Q is arbitrary if d(k) 
does not depend on kz- Since Imx° (Q,tf) is dominated 
by the contribution from the one dimensional bands, we 
consider the dynamical spin-susceptibility only for the 
one-dimensional bands in this paper. 

For the normal state with perfectly nested Fermi sur- 
face (^k' — "Ck) the imaginary part of the dynamical 
susceptibility is 



ImxL™ai(Q,^)-fA^(0) (1-2/(1 



(20) 



where f{uj/2) is the Fermi distribution function and the 
constant density of states N{0) is assumed. 

We study three possible cases for the triplet supercon- 
ductivity. In the first case (case A) the constant energy 
gap opens in the one-dimensional bands, while the line 
nodes are in the two-dimensional band. In the second 
case (case B) we assume that the order parameter in the 
one-dimensional bands depends only on as cos kz. In 
the last case (case C) the order parameter in the one- 
dimensional bands depends both on ky and kz and it is 
zero in the horizontal line nodes. 



k,k' ,(T,<7 



U' = ^ I 5ii{ sinfc^ sinfc^aj^^aL^_^a_k'-<T'aka 



sin fc„ sin fc!, 6,1.6^ ,, .^-k'-o-'^ko-'} 



y kcr — k— cr 

k' ki, 



+ 5i3{sinfc:rSin ycos^cos ^ 
^ (aLaLk-o-C-k'-CT'Ck'cr' + ft-.c.) 

7 • ^Z 

+ sm ky cos — sm — cos — 
^222 

^ (&L^-k-^c_k'-a'CkV' + h.c.)}^ , 



(24) 



and b! , are creation operators for the elec- 



where a^^ 

tron in one-dimensional bands and cj^^ is the creation 
operator for the two-dimensional band. In the above 
only terms relevant to the spin-triplet superconductivity 
with d(k) II z are included. The terms proportional to 
gii describe the attractive interaction between electrons 
with up and down spins in the one-dimensional bands 
which makes the triplet superconductivity with d vector 
parallel to the z axis. The 513 terms represent the pair 
hoppings between one-dimcnsional and two-dimensional 
Fermi surface. 

In this case Imx?j(Q,w) from the one-dimensional 
band is obtained as 



ImxL(Q,^) = 




—I — =, tanh ^ 



if w < 2A 
if w > 2A 

(25) 
if < 2A 



f A^(0) j^ -^ tanh ^ if > 2A 



(26) 
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FIG. 2: Imaginary part of the dynamical susceptibility nor- 
malized by (7r/2)7V(0) in case A (full gap on one- dimensional 
Fermi surface, eqs. 1211 and 1221 '). We take T = for the 
superconducting state and we take finite T for the normal 
state. (Since uj and T is scaled by A, we take T = O.IA for 
the normal state.) 



FIG. 3; Imaginary part of the dynamical susceptibility in case 
A (full gap on one-dimensional Fermi surface, eqs. 12111 and 
^) calculated in RPA. We take Ai = 0.001 eV, T = 0.0001 
eV, r = 0.0001 eV and U = 0.175 eV. 



where A = Ai sin kp is the energy gap on the Fermi sur- 
face in the one-dimensional bands. In Fig. [3 we plot 
imaginary part of the dynamical susceptibility normal- 
ized by (7r/2)iV(0) as a function of w at T = for the 
superconducting state and at a finite temperature for the 
normal state [T = O.IA). Since the coherence factor 

C'ir^(k, Q) is 1, Imx2^(Q,ti') diverges at w = 2A as the 
density of states in the s-wave superconductivity. 

In Fig. 121 we plot the dynamical susceptibility obtained 
in RPA with parameters, Ai = 0.001 eV, T = 0.0001 eV, 
F = 0.0001 eV and U = 0.175 eV for the one-dimensional 
band (Eq. lfT7|) and sinfe^ w 0.922). A peak appears at 
UJ = 2Ai sinfcF w 0.00184 eV only in ImXz^(Q,tj). 
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Next we study the case B. We assume horizontal line 
nodes on the one-dimcnsional Fermi surface as 

di2(k) = Ai sin /ca; cos (27) 
d22(k) = iAi sin fcy cos fc^ (28) 

In this case the Cooper pairs in the one-dimensional 
bands are formed between electrons on r and r + 
(0,0, ±c). Although this order parameter is not hkely 
to be realized in Sr2Ru04, we study this case to show 
the mechanism of resonance peak. In this case if we take 
Qz = or TT, line nodes are connected by Q as shown in 
Fig. ^ For = we obtain that 



FIG. 4: Nesting vector Q in the case B (horizontal line nodes 
on one-dimensional Fermi surface, eqs. 1271 and 12811 ). Thick 
solid and thin dashed lines show the side view of the Fermi 
surface on which order parameter is positive and negative, 
respectively. Solid arrows (Qz = 0) and dashed arrows [Qz = 
tt) are the nesting vectors. 
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FIG. 5: Imaginary part of the dynamical susceptibility in the 
case B (line nodes in one-dimensional Fermi surface eqs. (1271 
and 12811 ). We take the same parameters as in Fig. |5| 



ImxL(Q,^) 

N{Q)F (arcsin(^)|(2A)2)tanhi^ 



Imx1_(Q,c^: 
N{Q)E [. 



2A^ 



if w < 2A 
tanh ^ 
if > 2A 



(29) 



arcsm 



iV(0)i? (2^ 



if CJ < 2A 

tanh ^ 
if > 2A 



tanh 



4T 



(30) 



where F{ 



>|m) and E{4)\m) are incomplete elliptic inte- 
grals of the first and the second kinds, and K{m) and 
E{m) are complete elliptic integrals of the first and sec- 
ond kinds, respectively. As shown in the top panel of 
Fig. O ImXz2(Q, 1^) diverges at w = 2 A as the density of 
states for superconductivity with line nodes. For — tt, 
Imx^^ and Imx5|__ are exchanged. 
For = it/2 we get at T = 

ImxL(Q,^) = Imx" 


N{Q)F (arcsin((f 



K 



-(Q,^) 
if w < A 

if A < < \/2A (31) 



if > x/2A. 



In the RPA we get essentially the same result as shown 
in Fig. El 

For ^ 0, 7r/2, and tt, |di^(k)| Mi^(k')| at K = kp 
has two local maximum of 2Ai sin fci?| cos(Q2/2)| and 



FIG. 6: Imaginary part of the dynamical susceptibility in the 
case B (line nodes in one-dimensional Fermi surface eqs. (I27II 
and 12811 ) calculated in RPA. We take the same parameters 
as in Fig. El 
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FIG. 7: Imaginary part of the dynamical susceptibility in the 
case B (line nodes in one-dimensional Fermi surface eqs. 1271 
and 12811 ) calculated in RPA. We take the same parameters 
as in Fig. 



2AisinfcF|sin(Qj2)|. U Re{dl^(k)di^{k')) > at k^, 
where |(ii2(k)| -I- |diz(k')| becomes local maximum, then 

a peak appears in Imx-i (Q,a;), else a peak appears 

in ImXzziQ,^)- In Fig. [71 we plot the RPA result for 
Q = 7r/4. A peak is seen a.t u = 2AiC0S7r/8 w 
0.0017 in ImXzz(Q,w) and a small peak is seen a.t lu — 
2AiCOS7r/8 « 0.0007 in Imx+-(Q,w). 
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FIG. 8: Imaginary part of the dynamical susceptibility in the 
case C (line nodes on one-dimensional Fermi surface, eq. 1321 1 ). 
We take Ai = 0.001 eV, T = 0.0001 eV, F = 0.0001 eV. 



FIG. 9: Imaginary part of the dynamical susceptibility in the 
case C (line nodes on one-dimensional Fermi surface, eq. i'SZti ) 
calculated in RPA. We take the same parameters as in Fig.|3 



case C 



IV. COMPARISON WITH SPIN-SINGLET 
SUPERCONDUCTIVITY 



Finally we study the case C. Here we assume the order 
parameter in the one-dimensional bands to be 



di.(k) =d2.(k) 
"2 



Ai sin cos ■ 



y 



I cos — sm — 
2 2 



(32) 



This type of superconductivity is realized if the two- 
dimensional band is active and the one-dimensional 
bands are passivei^. In this case the F and Z points 
in Fig^ are not equivalent and the order parameter is 
zero at kz — ±7r, ±37r. We calculate Imxii(Q,ct;) nu- 
merically and plot Imxii(Q, cj) as a function of lo in Fig. 
IHl The small coherence peak exists at w « l.SAi. In 
this case Re((i*(k)(i2(k -)- q)) changes sign for k on the 

Fermi surface. The small coherence peak in (Qi^) 

for Qz = and that in Xzzi^,^) for Qz — 27r can 
be understood as follows. At — kp ^ O.GBtt and 
Q = (27r — 2kp, 27r — 2kp, Qz), the maximum value of 

|di.(k)| -f |di,(k')| 
— ^1 cos ^1^1 — COS kp cos ky 

+ I cos ^ k/l — cos kp co8{ky — 2kpy\ (33) 



is approximately 1.52Ai at ky « 1.627r and kz — 0, if 
= or 2tt. Since Re(dt^(k)di^(k')) > (< 0) at ky « 
1.627r and kz ^ for Qz ^0 {Qz = 2tt), linx+-(zz){Q,^) 
has a peak at w ~ 1.52Ai. 



Eq. Hll|) can be also applied for the spin-singlet 
case, when the order parameter is written as AQ^(k) — 
«fT^^A(k). As expected, we get the isotropic coherence 
factor 



/ , "aa'" 13/3''-' aa' 1313' 
aa'f3f)' 

eka'+Re(A*(k)A(k')) 



<Sy ( 1 ± 



(34) 



A peak appears in ImxM(q, w), only if Re(A*(k)A(k')) < 
0, which is the case for the d-wave paring with q w 
{■n,'K,qz) in high cupratesSl^LaSi^iSSii^. Note that 
the dynamical spin-susceptibility x°_(q, w) (eq. ((T5|) ') 
for the spin-triplet superconductivity is similar to 
the dynamical c/iar(7e-susceptibility in the spin-singlet 
superconductivitji^. 



V. CONCLUSION 

We study the dynamical spin-susceptibility for the 
spin-triplet superconductivity. The resonance peak ap- 
pears either in lTllXzz{^.,^) or Imx_| {^.,^) from the 

nesting part of the Fermi surface. We have shown that 
the existence of the line nodes on the quasi-one dimen- 
sional Fermi surface drastically changes the dynamical 
susceptibility, which can be observed by inelastic neutron 
scattering experiments. 

Coherence peak appears in Imx_| (Q,a;) if 

Re(d*(k)d^(k + Q)) > or in ImXz^(Q,w) if 
Re(d*(k)d^(k -I- Q)) < 0. When Re(d*(k)d^(k -I- Q)) 
changes sign in the nested part of the Fermi surface, 
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a coherence peak can appear in both Imxzz(Q,w) and 

Imxn (Q,^), but the divergence becomes weaker and 

easily smeared out. The position of the hne nodes on 
the Fermi surface, if hne nodes exist in the nested part 
of the Fermi surface, wiU be observed by scanning Qz in 
the inelastic neutron scattering. 

Recently, the resonance peak was searched in Sr2Ru04 
by inelastic neutron scattering^^iM, but no changes have 
been observed below Tc yet. We take the amplitude 
of the order parameter in the one-dimensional band as 



Ai = 0.001 eV in the RPA calculation. If the one- 
dimensional band is passive, the value may be smaller 
and the observation will be difficult. The gap structure, 
however, will be observed by experiments with better res- 
olutions using better single crystals. 

Recently, Mukuda et a^ reported a different depen- 
dence of the relaxation in Ruthenium and Oxide nuclear 
quadrupole resonance. This difference may be explained 

by the difference of xh (Qi'^) and Xzz{^.,^) in the spin 

triplet superconductivity. 
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